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ABSTRACT 

Using the Livermore X-ray Spectral Synthesizer, which calculates spectral models 
of highly charged ions based primarily on HULLAC atomic data, we investigate the 
temperature, density, and photoexcitation dependence of the /(II. 92 A)//(11.77 A) line 
ratio of Fe XXII. We find that this line ratio has a critical density ric ~ 5 x 10^^ cm~^, 
is approximately 0.3 at low densities and 1.5 at high densities, and is very insensitive 
to temperature and photoexcitation, so is a useful density diagnostic for sources like 
magnetic cataclysmic variables in which the plasma densities are high and the efficacy of 
the He-like ion density diagnostic is compromised by the presence of a bright ultraviolet 
continuum. Applying this diagnostic to the Chandra High Energy Transmission Grating 
spectrum of the intermediate polar EX Hya, we find that the electron density of its 
Te ~ 12 MK plasma is n^. = l.O^'jQg x 10^^ cm~'^, orders of magnitude greater than that 
typically observed in the Sun or other late-type stars. 

Subject headings: atomic processes — binaries: close — stars: individual (EX Hydrae) 
— stars: magnetic fields — X-rays: binaries 



1. Introduction 



Magnetic cataclysmic variables (mCVs; polars and intermediate polars) are bright X-ray 
sources because the highly supersonic [free-fall velocity = (2GMwd/-Rwd)^^^ ^ 3600 km s~^] 
material raining down on the white dwarf primary passes through a strong shock where most of 
its kinetic energy is converted into thermal energy and the plasma is heated to a temperature 
Tg = 3GM„d/^?77-H/8fci?„d ~ 200 MK. Because of the magnetic funneling of the mass lost by the 
secondary, the factor-of-four density jump across the accretion shock, and the settling nature of the 
post-shock flow (wherein the pressure is roughly constant and the density scales inversely with the 
temperature; Frank, King, &: Raine 1992), the hot plasma in mCVs is expected to be dense. For a 
mass-accretion rate M = 10^^ g s~^ (Lx = GM^^M /2R^^ « 3 x 10'^^ erg s~^) and a relative spot 
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size / = 0.1, the density of the flow immediately behind the shock is n = M/ATrfR^^fimii{vs/4) ^ 

The standard density diagnostic of high-temperature plasmas is the intensity ratio R = f/i of 
the forbidden to intercombination lines of He-like ions (Gabriel & Jordan 1969; Blumenthal, Drake, 
& Tucker 1972; Porquet et al. 2001), but this diagnostic is compromised in mCVs for two reasons. 
First, the critical density of this ratio increases with Z and hence temperature — the opposite of the 
trend in the accretion column — so the R line ratio is effective over only a very narrow slice of the 
accretion column. Second, in mCVs and other ultraviolet-bright stars, photoexcitation competes 
with collisional excitation to depopulate the upper level of the forbidden line, so the R line ratio 
can appear to be in the "high-density limit" if the radiation field is sufficiently strong at the 
appropriate wavelengths: for a plasma illuminated by a 30 kK blackbody, the R line ratios of all 
elements through Mg lie in the high-density limit, regardless of the density (Mauchc 2002). 

To circumvent the density/photocxcitation ambiguity of the R line ratio of Hc-Iikc ions, we un- 
dertook an investigation of potential density diagnostics of lines of Fe L-shclI ions that are observable 
in the high quality Chandra High Energy Transmission Grating (HETG) spectrum of the interme- 
diate polar EX Hya (Mauche 2000, 2002). Toward this end, Mauche, Liedahl, & Fournier (2001) 
investigated the density, temperature, and photoexcitation dependence of the 1(17.10 A)/7(17.05 A) 
line ratio of Fe XVII, and found that the anomalous ratio observed in EX Hya can be explained 
if the electron density of its Tc ~ 4 MK plasma is ric ^ 2 x 10^^ cm~^ or if the photoexcitation 
temperature Tbb ^ 55 kK. In this Letter, we investigate the density, temperature, and photoexci- 
tation dependence of the 1(11.92 A)//(11.77 A) line ratio of Fe XXII, which also is anomalous in 
EX Hya. We find that this line ratio is very insensitive to temperature and photoexcitation and 
that the observed ratio can be explained if the electron density of EX Hya's Tg ~ 12 MK plasma is 
ne ~ 1 X 10^^ cm~^. 

2. Model Spectra 

The density sensitivity of the X-ray lines of Fe XXII have been discussed previously by Doschek, 
Meekins, & Cowan (1973), Mason & Storey (1980), Doschek, Feldman, & Cowan (1981), Philhps 
et al. (1982), Fawcett et al. (1987), Philhps et al. (1996), and Wargehn et al. (1998) in studies of 
the X-ray spectra of solar and laboratory high-tcmperature plasmas. These diagnostics have not 
received more attention because of the relative weakness of Fe XXII emission lines (the Fe XXII 
ionization fraction peaks at just 22% and is greater than 1% over the relatively narrow temperature 
range Tc ~ 7 26 MK; Mazzotta et al. 1998), the relative weakness of n ^ 2 lines for n > 4, 
and the fact that the critical densities of these transitions are of order 10"*^^ cm^'\ larger than 
the densities typically observed in the Sun, tokamaks, and electron beam ion traps. We became 
interested in the density sensitivity of the 3^2 lines of Fe XXII because the critical density of 
this transition is comparable to the lower limit of the electron density derived from the Fe XVII 
7(17.10 A)//(17.05 A) line ratio in the Chandra HETG spectrum of EX Hya, and because the 
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Fe XXII 3^2 1(11.92 A)/7(11.77 A) line ratio is anomalous (the 11.92 A line is unusually strong) 

in the Chandra HETG spectrum of EX Hya compared to that of the bright late-type binaries 
HR 1099, Capella, and CrB (Ayres et al. 2001, Phillips et al. 2001, and Osten et al. 2003, 
respectively). 

Before wading into the details, we refer the reader to Figure 1, which illustrates the origin 
of the density dependence of the Fe XXII spectrum. At low densities (ng ^ 10^^ cm~^), the 
Fe XXII electron population is primarily in the 2s^2pi/2 ground state, and coUisional excitations 
are predominantly from the ground state into the 2s2p^ manifold (8 levels total) and the 2s^3d3/2 
level, both of which decay primarily to ground, producing lines in the extreme ultraviolet (EUV) 
(AA = 101-217 A) and at 11.77 A, respectively. However, these levels also decay to the 2s^2p3/2 
first-excited level with approximately 15% probability, producing lines in the EUV (AA = 114- 
349 A) and at 11.93 A; when that level decays to ground, it produces a line in the ultraviolet (UV) 
(A = 846 A). As the density increases {ue ^ 10^^ cm""^), electron population builds up in the 
2s^2p3/2 first-excited level because the Ml transition to ground is slow {A21 = 1.47 x 10^ s~^). At 
high densities, the 2s'^2p^/2 first-excited level is fed primarily by radiative decays from the 2s2p^ 
manifold. CoUisional excitations out of the first-excited level are primarily into the 25^3^5/2 level, 
which decays primarily back to the first-excited level, producing a line at 11.92 A. Consequently, 
the Fe XXII 11.92 A line is relatively strong in the X-ray spectra of high-density plasmas. 

To calculate quantitative models of the X-ray spectrum of Fe XXII, we employed the Livermore 
X-ray Spectral Synthesizer (LXSS), a suite of IDL codes that calculates spectral models of highly 
charged ions based on Hebrew University/Lawrence Livermore Atomic Code (HULLAC) atomic 
data. HULLAC calculates atomic wavefunctions, level energies, radiative transition rates Ay^, and 
oscillator strengths /y according to the fully relativistic, multiconfiguration, parametric potential 
method (Klapisch 1971; Klapisch et al. 1977). Electron impact excitation rate coefficients arc 
computed quasi-relativistically in the distorted wave approximation (Bar-Shalom, Klapisch, & Orcg 
1988) assuming a Maxwellian velocity distribution. Our Fe XXII model includes electron impact 
excitation rate coefficients and radiative transition rates for El, E2, Ml, and M2 decays for levels 
with principal quantum number n < 5 and azimuthal quantum number Z < 4 for a total of 228 
levels. Using these data, LXSS calculates the level populations for a given temperature and density 
assuming collisional-radiative equilibrium; the line intensities are then simply the product of the 
level populations and the radiative decay rates. 

In a preliminary report, we (Mauche, Liedahl, Sz Fournier 2003) used LXSS to calculate the 
X-ray spectrum of Fe XXII for electron densities Uc = 10^^ to 10^^ cm^^ and electron temperatures 
Te = 6.3, 12.8, and 25.5 MK (~ ^, 1, and 2 times the temperature at which the Fe XXII ionization 
fraction peaks; Mazzotta et al. 1998). At the resolution of the Chandra HETG, the 2s'^3d2/2~ 
2s'^2pi/2 line at 11.77 A is blended with the {very much weaker) 2s'^3d^/2~^^'^^Pi/2 line at 11.76 A, 
and the 2s^ 3^5/2-25^2^3/2 line at 11.92 A is blended with the (much weaker) 2.5^3^3/2^25^2^3/2 line 
at 11.93 A, so we report as 7(11.92 A)//(11.77 A) the ratio of the sum of these blends. Figure 3 of 
Mauche, Liedahl, Sz Fournier (2003) shows that this line ratio has a critical density Hc ~ 10^^ cm~^. 
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is approximately 0.3 at low densities and 1.4 at high densities, and has only a weak temperature 
dependence. 

2.1. Modifications to LXSS: CoUisional Excitation 

Motivated by the detailed study of Fournier et al. (2001) of the density-sensitive Fe lines in 
the EUV spectrum of the high-density Frascati Tokamak Upgrade plasma, for this Letter we in- 
vestigated the effect on the Fe XXII 1(11.92 A)//(11.77 A) line ratio of two modifications to the 
collisional excitation data used in LXSS. First, we investigated the effect of replacing, for all transi- 
tions between and among the 2s^2p and 2s2p^ levels of Fe XXII, the electron impact excitation rate 
coefficients computed with HULLAC with those of Zhang & Pradhan (1997), computed with the 
relativistic R-matrix method. At Tc = 12.8 MK, the Zhang & Pradhan electron impact excitation 
rate coefficients are larger than those of HULLAC by 1.95 for the 2s^2px/2^2s^ 2^3/2 transition, 
0.92-1.70 (average 1.27) for the eight 2s'^2pi/2-2s2p'^ transitions, 1.05-3.23 (average 1.55) for the 
eight 2s^2p3/2-2s2p^ transitions, and 0.45-8.99 (average 2.20) for transitions between levels within 
the 2s2p^ manifold (but only 1.00-1.02 for the four 2s'^2p-2s'^3d transitions). Despite the significant 
differences between these excitation rate coefficients, we find that for Tc = 12.8 MK the Fe XXII 
7(11.92 A)//(11.77 A) line ratio increases by 7% at Uc = 10^^ cm^^, but the low- and high-density 
limits are unaffected. Second, we investigated the effect of adding proton excitations for transi- 
tions among the levels of the 2s^2p and 2s2p^ configurations. This was accomplished by fitting the 
proton impact excitation rate coefficients of Foster, Keenan, & Reid (1997) to an expression of the 
form Q{Te) = 10"^^ Te~^^^ Yl^=i A exp(— fij/Te), where Q{Te) is the proton impact excitation rate 
coefficient in units of cm^ s~^, and we have assumed that the proton and electron temperatures are 
equal. Adding proton excitations to the LXSS population kinetics calculation under the assumption 
that the proton number density Hp = 0.85 ne (as is appropriate for a high-tcmpcrature plasma with 
solar abundances), we find that for Tg = 12.8 MK the Fe XXII /(11.92 A)//(11.77 A) line ratio 
increases by 19% at = 10^^ cm~^, and the high-density limit increases by 7% to approximately 
1.5. 

With both of these modifications to LXSS, we recalculated the X-ray spectrum of Fe XXII for 
electron densities Ue = 10^^ to 10^^ cm~^ and electron temperatures Te = 6.3, 12.8, and 25.5 MK. 

The upper panel of Figure 2 shows that the model Fe XXII 7(11.92 A)//(11.77 A) line ratio has a 
critical density nc ~ 5 x 10^^ cm~^, is approximately 0.3 at low densities and 1.5 at high densities, 
and has only a very weak temperature dependence. 

2.2. Modifications to LXSS: Photoexcitation 

To investigate the effects of photoexcitation on the Fe XXII 7(11.92 A)/7(11.77 A) line ratio, 
we added to the LXSS population kinetics calculation the photoexcitation rates (ne^ /mec)fijF^{T), 
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where Fi,{T) is the spectrum of the photoexcitation radiation field. For simphcity, we assume that 
Fy[T) = [Att /hu)By{T\^\^) (i.e., the radiation field is that of a blackbody of temperature T^b) and 
the dilution factor of the radiation field is equal to ^ (i.e., the X-ray emitting plasma is in close 
proximity to the source of the photoexcitation continuum). With these assumptions, we used LXSS 
to calculate the X-ray spectrum of Fe XXII for an electron temperature = 12.8 MK, electron 
densities ng = 10^^ to 10^^ cm~^, and photoexcitation temperatures = 20, 40, 60, 80, and 100 
kK. The lower panel of Figure 2 shows that below Tbb = 60 kK the Fe XXII /(11.92 A)/7(11.77 A) 
line ratio has essentially no photoexcitation dependence, but as the photoexcitation temperature 
increases the line ratio increases at low electron densities and the critical density shifts to higher 
densities, while the high-density limit is unaffected. Compared to the Fe XVII /(17.10 A)//(17.05 A) 
line ratio, the Fe XXII 7(11.92 A)//(11.77 A) hne ratio is insensitive to photoexcitation because 
the 2s^2pi/2~2s^2p3/2 transition in the UV is not optically allowed (/12 = 3.15 x 10~^) and the 
2s^2pi/2~2s2p^ transitions lie at such short wavelengths in the EUV. 

3. Observations 

The Chandra HETG/ACIS-S observation of EX Hya was performed between 2000 May 18 
9*^41™ and May 19 2'^^54™ UT for a total exposure of 59 ks. Extraction of the grating spectra 
and calculation of the effective area files was accomplished with the CIAO 2.1 suite of software 
using the reprocessed data products and new calibration data files (version R4CU5UPD8) for 
sequence 300041. Various aspects of this spectrum have been discussed by Mauche (2000, 2002) 
and Mauche, Liedahl, &; Fournier (2001); here we discuss only a tiny portion (AA 11.60-12.05 A) 
of the medium-energy grating (MEG) spectrum as it bears on the relative strengths of the Fe XXII 
n = 3 ^ 2 emission lines. This spectrum is shown in Figure 3, where the ordinate is in units of 
counts s~^ A~^, the wavelength bin width AA = 0.005 A, and we have combined ± first order 
counts. This portion of the spectrum contains the 11.77 and 11.92 A emission lines of Fe XXII as 
well as the 11.74 A emission line of Fe XXIII. 

To determine the flux in these emission lines, we fit the MEG spectrum over the AA = 11.60- 
12.05 A wavelength interval with a model consisting of a linear background (to account for the ther- 
mal bremsstrahlung continuum obvious in the broadband spectrum) and three Gaussians. Assum- 
ing a common Gaussian width, the fitting function /(A; a) = oi + a2(A — Aq) + X^^=o '^5+2i exp[— (A — 
'3^4+2i)^/2a3]/\/27ra3. To perform the fit, we separately accounted for it first orders and used 
AA = 0.005 A bins, so there were 180 data points and 171 degrees of freedom (dof). The 
resulting fit, combining zb first orders, is shown by the thick histogram in Figure 3. It gives 
X^/dof = 104/171 = 0.60 and fit parameters as follows: Gaussian width 03 = 1.41q g mA and line 
wavelengths and fluxes 7(11.740 A) = 7.74 (±0.80) x 10-^ 7(11.776 A) = 4.27 (±0.67) x 10-^ and 
7(11.923 A) = 4.52 (±0.68) x 10"^ photons cm s , where the errors on the wavelengths are ± 
1-2 mA and all errors are 1 a (68% confidence) for one interesting degree of freedom. Formally (e.g., 
Bevington 1960), the 7(11.92 A)/7(11.78 A) line ratio equals 1.06±0.23 (68% confidence), but this 
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propagation of errors ignores a subtle interplay between the error bars of the numerator and de- 
nominator of this ratio. To generate more rigorous values for the error bars of this ratio, we used a 
Monte Carlo technique with 10® realizations to determine that the observed 7(11.92 A)//(11.78 A) 
line ratio equals 1.06 and the errors are (68% confidence), ^032 (90% confidence), and I047 

(99% confidence). This best-fit line ratio and these error envelopes are shown superposed on both 
panels of Figure 2. 



4. Analysis 

Figure 2 shows that the combination of our LXSS models and Chandra HETG spectrum places 
significant constraints on the density of the Tq 12 MK plasma in EX Hya. Specifically, allowing for 
a factor-of-two uncertainty in the electron temperature, the upper panel of Figure 2 shows that the 
electron density rie = I.O^q'j x 10^^ cm^"^ at 68% confidence, Uc > 3.2 x 10^'^ cm^^ at 90% confidence, 
and Tie > 1.7 X 10^^ cm~^ at 99% confidence. Conversely, the lower panel of Figure 2 shows that 
the electron density can be arbitrarily low if the photoexcitation temperature Tbb ^ 100 kK, but 
that alternative is excluded by UV observations of EX Hya. First, assuming that the distance to 
EX Hya is 100 pc, that the radius of its white dwarf is i?wd = 10^ cm, that its 1010 A flux density 
equals 2.5 x 10^^^ erg cm^^ s~^ A^^ (Mauchc 1999), and that 100% of this emission is due to a 
100 kK hotspot (i.e., assuming that the white dwarf, accretion disk, and accretion curtains do not 
radiate in the UV), tlic fractional emitting area of the spot is approximately 0.5%, far smaller than 
is inferred from optical and UV light curves of mCVs. Second, fits to EX Hya's UV spectrum yield 
effective temperatures T^s ~ 25 kK (Greeley et al. 1997; Eisenbart et al. 2002), far lower than that 
required to produce the anomalous Fe XXII 7(11.92 A)/7(11.77 A) line ratio via photoexcitation. 



5. Discussion 



To date, four different spectroscopic diagnostics have provided evidence of high densities in 
the X-ray-emitting plasma of the intermediate polar EX Hya. Hurwitz et al. (1997) used the line 
ratio of the Fe XX/Fe XXIII 133 A blend to the Fe XXI 129 A line observed in the 1994 Extreme 
Ultraviolet Explorer spectrum of EX Hya to infer rig ^ 10^^ cm~^ under the assumption that 
the plasma temperature = 10 MK. Using the Chandra HETG spectrum of EX Hya, Mauche 
(2002) showed that the He-like ion R line ratios of O, Ne, Mg, Si, and S are all in their high- 
density limit; Mauche, Liedahl, & Fournier (2001) used the Fe XVII 7(17.10 A)/7(17.05 A) line 
ratio to infer ric ^ 2 x 10^^ cm~^ for its Tg ~ 4 MK plasma; and we here used the Fe XXII 
7(11.92 A)/7(11.77 A) line ratio to infer rig « 1 x lO^'' cm'^ for its Tg « 12 MK plasma. Of 
these diagnostics, the Fe XXII line ratio is the most reliable because it has the highest critical 
density and is the least sensitive to temperature and photoexcitation. We conclude first that the 
Chandra HETG spectrum of EX Hya requires plasma densities that are orders of magnitude greater 
than those typically observed in the Sun or other late-type stars, and second that the Fe XXII 
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7(11.92 A)/7(11.77 A) line ratio [like the Fe XVII 1(17.10 A)/7(17.05 A) line ratio] is a useful 
density diagnostic for sources like mCVs in which the plasma densities are high and the efficacy 
of the He-likc ion density diagnostic is compromised by the presence of a bright UV continuum. 
Finally, we note that our density determinations for EX Hya are consistent with those expected 
for plasma in the accretion column, where rie 9S T~^. Additional density diagnostics are needed to 
begin to map the temperature/density profile of the accretion column of this and other mCVs. 
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under contract NAS8-39073. This work was performed under the auspices of the U.S. Department 
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Fig. 1. — Schematic diagram of the dominant colhsional and radiative processes of Fe XXII for 
electron densities rie = 10^^ (a) and 10-*^^ cm~^ (6). Collisional excitations are indicated by upward- 
pointing straight hnes; radiative decays are indicated by downward-pointing wavy hnes and are 
labeled with the wavelength of the transition in Angstroms. 
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Fig. 2.— LXSS model Fe XXII 7(11.92 A)//(11.77 A) line ratio as a function of electron density 
for (a) electron temperatures = 6.3, 12.8, and 25.5 MK {dotted, solid, and dashed curves, 
respectively) and (6) electron temperature = 12.8 MK and photoexcitation temperatures Ty^ = 0, 
60, 80, and 100 kK (solid, dotted, short dashed, and long dashed curves, respectively). Horizontal 
line and the dark, medium, and light shaded stripes indicate respectively the value and 68%, 90%, 
and 99% confidence error envelopes of the line ratio measured in EX Hya. 
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Fig. 3. — Detail of the Chandra MEG spectrum of EX Hya in the neighborhood of the Fe XXII 3^2 
Unes. Data are shown by the filled circles with error bars and the net and separate contributions 
to the model fit are shown by the histograms. Data combines it first orders and are binned to 
0.005 A. 



